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Abstract 

Using Monte Carlo simulation we investigated time of flight current transients 
predicted by the dipolar glass model for a random spatial distribution of hopping 
centers. Behavior of the carrier drift mobility was studied at room temperature 
over a broad range of electric field and sample thickness. A flat plateau fol- 
lowed by j oc t~ 2 current decay is the most common feature of the simulated 
transients. Poole-Frenkel mobility field dependence was confirmed over 5 to 
200 V/fim as well as its independence of the sample thickness. Universality of 
transients with respect to both field and sample thickness has been observed. 
A simple phenomenological model to describe simulated current transients has 
been proposed. Simulation results agree well with the reported Poole-Frenkel 
slope and shape of the transients for a prototype molecularly doped polymer. 
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1. Introduction 

Dipolar glass (DG) model 1-5] was developed in the late 1990s in response 
to the urgent need to explain the ubiquitous Poole-Frenkel (PF) mobility field 
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dependence observed in amorphous polymers (polyvinylcarbazole [6(] , polysilenes 
[7] and others), molecularly doped polymers |§4llj]. an d low molecular weight 
organic glasses [l2[ in a broad field range. This model overcomes the limitations 
of the Gaussian Disorder Model (GDM) [l3j], which successfully explained the 
mobility temperature dependence, but failed to reproduce the mobility field 
dependence for small and moderate electric fields. The DG model is a natural 
evolution of the earlier approach of Borsenberger and Bassler [Til HBj . who 
suggested dipolar energetic disorder as an explanation of strong polarity effect on 
the hopping charge mobility in amorphous organic materials. They considered 
the dipolar as well van der Waals energetic disorder, but did not introduce 
them consistently on a microscopic level. It turned out that strong spatial 
correlations in the random energy landscape, intrinsic to dipolar disorder, move 
transport properties of polar materials far away from those of the GDM. In 
particular, spatial correlations naturally provide a strong physical foundation 
for the development of the PF mobility field dependence [2, 3]. 

Previous Monte Carlo (MC) simulations [16| have shown that the time of 
flight (TOF) currents, predicted by the DG model, provide transients having 
typical features of the experimental ones: an initial short spike, reflecting a non- 
equilibrium stage of the carrier transport, followed by a flat plateau signaling 
partial equilibration of carriers, and then an anomalously long tail described by 
a power law j oc (/3 w 2.0 — 2.5). All earlier models failed to predict these 
broad post flight current tails in combination with almost flat plateau for the 
case of strong disorder. 

The aim of the present paper is to investigate in detail the TOF currents, 
predicted by the DG model for a prototype molecularly doped polymer 30% 
DEH:PC (polycarbonate doped with 30 wt.% of aromatic hydrazone DEH). 



2. Basics 

We consider charge transport in a cubic lattice of randomly oriented static 
dipoles with the dipole moment p, where c is the fraction of sites serving as 
hopping centers (transport sites) . The random energy of a particular transport 
site is Ui = &p(ri), where tp(fi) is the electrostatic potential, created by all other 
dipoles. If c is not too low, then the density of states has a Gaussian form and 
correlation function C{r) = {U(r)U(0)} decays as 1/r [l|, 0, [l7j]. If electric field 
F is applied, then the random energy has an additional term —eFri. The DG 
model is also known as a Correlated Disorder Model (CDM) Q, but the present 
name is better suited for description of the true nature of the model; for example, 
it explicitly assumes a particular spatial decay of the correlation function C(r) oc 
1/r, while the name "CDM" is too ambiguous in this respect. Indeed, the model 
of the quadrupolar glass, where randomly oriented quadrupoles fill the lattice 
instead of dipoles, provides a correlated energy landscape as well, so it again 
can be dubbed a Correlated Disorder Model, but the correlation function in the 
quadrupolar case is different, CJr) oc 1/r 3 [3 [3, and, hence, the transport 
properties are different as well [20j. 
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The Miller-Abrahams hopping rate was used for the simulations, where the 
rate of transition from site i to site j is given by 



fexpf-^r^V U i -U i >0 
Pi^j = l> exp(-27rij) < V kl J (1) 
[l, Uj-Ui<0 

here vq is the prefactor frequency, r,-j = \fj — f^|, and 7 is a wave function decay 
parameter for transport sites. 

Monte Carlo simulations refer to the traditional TOF experiment when a 
sheet of carriers (holes in the case of DEH:PC polymer) is instantly produced 
at the generating electrode at x = and then drifts towards the collecting 
electrode located at x — L (hence, L is a thickness of the transport layer). 
The output of the simulation is the mean velocity v (t) of a carrier; the major 
experimental observable, current density j(t), is proportional to v(t). Another 
important parameter is the mean carrier velocity (v) for a full transfer from 
the generating to collecting electrode. This velocity is equal to (L/td), where 
id is the time for a carrier to drift from x = to x — L. Other details of the 
simulation can be found elsewhere (3. 13|. 



In the simulation we tried to model the polycarbonate doped with 30 wt.% of 
aromatic hydrazone DEH. DG model parameters (borrowed from the analysis 
of experimental data in Ref. [2l|) are as follows: the mis dipolar disorder a 
is 0.13 eV and contains no additional contributions (e.g., van der Waals), the 
lattice constant a — 0.77 nm, 27a = 11.8, kT = 0.0252 eV (room temperature), 
and c = 0.3. 

In all figures (apart from Fig. I10[) we use a dimensionless time t = vq exp (—27a) t' , 
where t' is the physical time. Frequency v$ was taken to be 1 x 10 16 s _1 to pro- 
vide matching of the simulated carrier velocity at 5 V//mi with the experimental 
one from Ref. [2l| . 



3. Simulation results 

Fig. Q] shows TOF transients for L = 15.4 /mi (20 000 lattice planes) for 
fields in the range 5 — 200 V///m. In all cases the transport occurs in a quasi- 
stationary regime with well-defined transient plateaus. Log-log representation 
shows that the post flight current decay follows a power law j oc t~@ 2 rather 
than an exponential one, usually expected for fully equilibrated transport with 
well-defined plateau of the transient. Time of flight tt r is determined as the 
intersection of the tangents to the initial (slope /3i) and post flight (slope P2) 
parts of the curve as shown by broken straight lines in the figure. While /3i is 
close to zero, parameter /?2 varies between 2.1 and 2.7. For fields exceeding 36 
V//zm the initial current spike begins to appear indicating a non-equilibrium 
stage of the transport process. 

Current transients can be presented in double linear coordinates to give tra- 
ditional transit times to or ti/ 2 (shown in Fig. [2|). Parameter W = (ti/2 — to) A1/2 
is approximately constant in the range 5 — 200 V//im, fluctuating around 0.45. 
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Figure 1. Simulated TOF transients for the electric field equal to 5 (f ), 10 (2), 
36 (3), 80 (4), 110 (5), 150 (6), and 200 V/^m (7), correspondingly; L = 20 000 
lattice planes. Time of flight t tr is shown by the arrow. 

Transit times (t tr , to, and £1/2) for various fields are summarized in Table [TJ 
One can see that U r is bounded by t-yji (upper bound) and to (within uncertain- 
ties in ttr due to errors in plotting the asymptote to the tail of the transient). 
Fig. [3]shows field dependence of the relative drift mobility [i = L/Ft, calculated 
in four possible ways, using the corresponding values of tt r , to, t\i<i, as well as 
MC drift time L/ (v) (shown is the ratio of the mobility to the corresponding 
mobility at 5 V//zm). For low and moderate fields all four curves are rather 
similar, they generally follow PF dependence 

Mo <x exp (SF 1 / 2 ) (2) 

where S is the PF slope. At high fields the calculated mobility field depen- 
dence begins to deviate from $2$ , reflecting the particular property of the Miller- 
Abrahams hopping rate [3j. 

Thus, three main results of the previous simulations (fast establishment of 
the quasi-stationary transport regime as evidenced by the flat plateau, the power 
law decay of the tail of the transients, and PF mobility field dependence) have 
been confirmed. 

A new observation concerns the carrier equilibration time t eq , which is easily 
seen in Fig. [T]at large electric fields but seems rather difficult to quantify as the 
non-equilibrium stage merges smoothly into the equilibrium one. To clarify this 
situation we present Fig. |4]- Fig. [6] showing TOF transients at extreme fields 
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Figure 2. Transient 1 from Fig. Q](F = 5 V//xm, L = 20 000 lattice planes) in 
linear coordinates with to and ti/ 2 indicated by arrows. 

(5 and 200 V//im, respectively) for varying thickness L. Equilibration time t eq 
can be estimated as the time of an intersection of the tangents to the initial 
non-equilibrium part of the transient and the plateau. 

The use of this procedure gives the equilibration time 1.4 x 10 6 for the highest 
field of 200 V//im, while for the smallest field (5 V/fim) it is much larger and 
equal to 1.2 x 10 8 . The observed variation of t eq is even larger than that of the 
carrier mobility: the former rises by a factor of 80, while the latter increases by 
a factor of 20. 

Fig. [l] Fig. |4j and Fig. [5] demonstrate features, usually attributed to 
the current universality traditionally observed in the continuous time random 
walk theory of Scher and Montroll [221 ] or multiple trapping formalism with the 
exponential trap distribution [23[ • Indeed, the pre- flight part of the TOF curves 
for the DG model is a flat plateau with /3i ss 0, while the post flight decay again 
follows the power law dependence with the exponent 02 ~ 2.0. These values 
agree well with the Scher-Montroll theory, where fi\ = 1 — a and 02 — 1 + a with 
a < 1. Additionally, the transients demonstrate universality with respect to the 
field and thickness variation. The only contradicting factor is the existence of 
the equilibration time, which does not scale properly with the time of flight. 
For this reason we expect a mild violation of universality with respect to F for 
high fields, as corroborated by Fig. [Jj 

Simulation data allows one to study the thickness dependence of the char- 
acteristic times of the transport process (Fig. [4]- Fig. [6j. Obviously, t eq does 
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Figure 3. Mobility field dependence (data points refer to the mobility defined 
by the different methods: using to (□), t 1 j 2 (■), t tr (°), and (v) (•), corre- 
spondingly (/15 is the corresponding mobility for F = 5 V//zm). Solid lines 
serve as a guide to the eye, and the broken line corresponds to the PF slope 
0.32 (/xm/V) 1 / 2 . The lowest curve (A) shows the GDM mobility, calculated 
using (v). 



not depend on L. Times of flight found on log- log as well as lin-lin plots exhibit 
thickness dependence very close to a linear one as could be seen from Fig. [8] 
and Table [5] (for F = 200 V//im and L = 1000 transient is too dispersive to 
calculate to and t 1 / 2 ). 

Until now we considered the classical TOF geometry where an instantaneous 
generation of the thin sheet of carriers takes place near the generating electrode. 
Recently two more geometries have been studied experimentally, namely, an 
instantaneous uniform generation of carriers in the bulk of the layer (TOF-2), 
or similar bulk carrier generation with a regulated depth of the generation zone 
(TOF-la) 24]. Fig. [9] compares all three TOF variants for an equal total 
number of generated carriers (as always, in a small signal regime). As expected, 
if the TOF transient demonstrates a flat plateau, then the corresponding TOF-2 
curve demonstrates an almost linear decay in the time range where the TOF 
current remains constant. Both transients decay similarly in the post-flight 
region, but the TOF current is always greater than the TOF-2 one, their ratio 
being practically constant (in our case it is equal to 3.2) as the log-log plot 
shows (Fig. [9l inset). 

The TOF-la transient, simulated for the thickness of the generation layer 
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Table 1. Field dependence of transit times for L = 20 000 lattice planes. 



F, V//im 


Time of flight 


io-\ r 


10~% 


10- 8 ii/2 


5 


346 


197 


490 


10 


144 


100 


187 


20 


41.1 


34.1 


62.7 


36 


15.6 


11.0 


22 


55 


6.78 


4.71 


9.1 


80 


3.3 


2.0 


4.2 


110 


2.2 


1.1 


2.4 


150 


1.1 


0.46 


1.1 


200 


0.5 


0.34 


0.66 



Table 2. Thickness dependence of transit times. 



F, V/^m 


L, lattice planes 


Time of flight 


10~V 


io-% 


10-%/2 


200 


1000 


0.0228 






50 000 


1.26 


0.925 


1.68 


100 000 


2.74 


1.96 


3.40 


200 000 


5.26 


4.04 


6.76 


5 


500 


9.05 


5.00 


12.10 


5 000 


85.1 


47.1 


119.0 


20 000 


346 


197 


490 



equal to 0.5L, demonstrates typical features of quasi-equilibrium transport. The 
plateau becomes more tilted and its length is approximately twice shorter (tak- 
ing into account the accuracy of the MC simulation), as expected in this case 



4. Discussion 

The most prominent features of the DG model are the PF field dependence 
of the carrier mobility, the particular shape of transients demonstrating the 
combination of a flat plateau with the slow post-flight current decay that follows 
a power law with the exponent 02 close to 2.0, and current universality with 
respect to F and L. The first feature was thoroughly discussed and explained 
earlier Here we would like to discuss the specific behavior of transients. 

It is worth noting that universality of transients has alread y b een reported in 



experimental papers for non-dispersive charge transport 26j-[29j. 

Evidently, in our case development of a well-defined plateau cannot be de- 
scribed by the usual approach where the long time behavior is governed by the 
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Figure 4. Simulated TOF transients for L = 500 (1), 5 000 (2), and 20 000 (3) 
lattice planes, correspondingly. Electric field is 5 V//j,m. 



diffusion equation (such as the multiple-trapping model in the case of a fast de- 



caying trap distribution |30|,|3l|) and current decay follows the exponential time 



dependence. An alternative explanation was suggested by Nikitenko et al. [321 ] . 
who calculated an effective diffusion coefficient for the GDM and found that it 
slowly increases with time on a scale much longer that the typical equilibration 
time of the mean carrier velocity. Detailed analysis of the long time behavior 
of transients was not performed but the particular case, shown in Figure 3 of 



Ref. [32j , demonstrates approximate power law decay of the transient with the 



exponent $<i close to 3. Applicability of this approach to the DG model is not 
evident, especially keeping in mind a well-established difference in the mobility 
field dependence in the GDM and DG model 0, [Hj]. In addition, MC sim- 
ulations indicate that in the GDM for the case of strong disorder a/kT ~ 5 
transients are much more dispersive: they do not develop a flat plateau and 
demonstrate a superlinear dependence of the transit time on the sample thick- 
ness 33, HJ. For these reasons we suggest a following simple phenomenological 
description of the charge transport. 

Formally, transport behavior found in our simulation can be described by 
the "quasi-ballistic" model (qBM) where a distribution of drift carrier velocities 
takes place in a trap free organic solid. Indeed, let f(v) be the density distribu- 
tion of equilibrated carrier velocities. If a carrier moves with some velocity v, 
then at time t — L/v it reaches the collecting electrode and does not contribute 
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Figure 5. Simulated TOF transients for L equal to 1 000 (1), 50 000 (2), 100 000 
(3), and 200 000 (4) lattice planes, correspondingly. Electric field is 200 V//nm. 



to the total current anymore. Hence, in this model the current transient is 

oo 

j(t) =en I dvvf(v)6 (L - vt) (3) 



o 



where 9 (x) is a unit step function, and n is the initial density of carriers. Cal- 
culating the derivative, we obtain 



J t =- en J dvv 2 f(v)5 (L - vt) = -en^-f (L/t) (4) 



or 

t(v) = ^ -r 



(5) 



t=L/v 



In general, if j(t) oc t~P 2 for t — > oo, then f(v) oc v^ 2 ~ 2 for v — > 0. Accordingly, 
for ft = 2 we have f(v) rts const for v — > 0. 

In this model TOF transients are universal and all transit times are strictly 
proportional to L. Assuming an abrupt decay of f(v) at v = v max we auto- 
matically obtain a flat plateau for t < L/v max . Also, TOF-2 transient should 
have a linear decay extending for a full length of the TOF plateau, while the 
TOF- la transient should have a shortened plateau in comparison to the original 



9 




Figure 6. Simulated TOF transients in linear representation for L equal to 
50 000 (1), 100 000 (2), and 200 000 (3) lattice planes, correspondingly. These 
transients are the same as curves 2 — 4 in Fig. [5j Electric field is 200 V/fim. 



TOF one. Thus, all dynamic features predicted by the DG model for all three 
variants of the time of flight technique are reproduced by the model. 

Limitations of the qBM are obvious: it cannot reproduce the non-equilibrium 
stage of the carrier transport, and the PF mobility field dependence has to 
be explicitly embedded into the model. Yet the real problem is to derive the 
phenomcnological ballistic model directly from the basic hopping transport of 
the DG model. 

At last we would like to discuss implications of the MC simulations regarding 
the prototype molecularly doped polymer 30% DEH:PC thoroughly investigated 
in literature [ill, 0, lipf . At room temperature, it has been reported that PF 
slope is equal to 0.39 |2l|j. 0.42, and 0.46 (/xm/V) 1 / 2 [Hj]. The PF slope for 
the straight line in Fig. [3] is slightly smaller and equals 0.32 (/um/V) 1 / 2 . 

Also, the simulated TOF transients are very similar to experimental ones 
[HI, they demonstrate the same flat plateau and slow power law decay (Fig. 
[TO]) . Systematic deviation at small t is probably connected with electrode effects 
(roughness of the electrode surface, presence of impurities, variation of the local 
structure of the material at the vicinity of the electrode). 
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Figure 7. Normalized transients for the electric field equal to 5 (1), 32 (2), 
110 (3), and 200 V//zm (4), correspondingly; L = 20 000 lattice planes. Curves 
3 and 4 are almost indiscernible in this presentation. Inset shows the analogous 
plot with ttr used as normalization time. Worse visual universality certainly 
follows from the much less reliable determination of t tT from log-log plots. 

5. Conclusions 

We performed Monte Carlo simulation for the DC model with a partially 
filled lattice where transport molecules occupy randomly 30% of all sites. Anal- 
ysis of simulated TOF transients confirmed the model's unique properties such 
as observation of Poole- Frcnkel mobility field dependence in a broad range from 
5 to 200 V/^m, development of a flat plateau {fix w 0) followed by the slow 
power law current decay {(3 2 ~ 2.0) and transient universality with respect to 
both field and thickness. 

It has been shown that due to this universality one may use either lin-lin 
(using to as well as tx/2 variant) or log- log plots to investigate the field or 
thickness dependence of the carrier mobility. 

Detailed examination of the TOF transients at strong fields reveals the exis- 
tence of the initial non-equilibrium relaxation of charge carriers. Equilibration 
time is typically short in comparison to the drift time, does not depend on the 
sample thickness and exhibits strong field dependence, even more pronounced 
than the corresponding field dependence of the mobility. 

Temporal behavior of the TOF currents (including TOF-2 and TOF-la) is 
well reproduced by a simple "quasi-ballistic" model which assumes a statistic 
(and time-independent) distribution f(v) of carrier drift velocities. If /(0) = 
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Figure 8. Thickness dependence of time of flight for t tr (B), to (°), an d ii/2 
(•), correspondingly; t(500) is the corresponding time for L = 500a and straight 
line represents the linear dependence t(L) = Li(500)/500. Electric field is 5 
V/fim. 

const > and f(v) drops abruptly at some velocity, then this model fairly well 
describes most prominent features of the simulated transients, the flat plateau 
followed by the power law decay j oc i~ 2 . 

An attempt to fit DG simulation data to a prototype molecularly doped 
polymer was successful regarding the Poole-Frenkel slope (0.32 versus 0.39 
(/um/V) 1 / 2 ), if measured at moderate electric fields, and shape of the transients. 
To the best of our knowledge, at the moment there is no alternative transport 
model capable of reproducing both the experimental mobility field dependence 
and shape of the transients for some particular organic material in the direct 
MC simulation using transport parameters (a, 70, and c), extracted from the 
TOF experiments for the same material. 
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Figure 9. Simulated TOF (1), TOF-la (2) and TOF-2 (3) transients in linear 
coordinates. L = 20 000 lattice planes, F = 5 V///m. Dotted straight line has 
been drawn to emphasize the linear decay of transient 3 in the time interval 
corresponding to the plateau of transient 1; tTOF-iaATOF ~ 0.63. Inset shows 
TOF (1) and TOF-2 (3) transients in log-log plot. 
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